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In the Minimal Supersymmetric Standard Model, the effective b quark Yukawa coupling to the lightest neutral 
Higgs boson is enhanced. Therefore, the associated production of the lightest Higgs boson with a b quark is an 
important discovery channel. We consider the SUSY QCD contributions from squarks and gluinos and discuss 
the decoupling properties of these effects. A comparision of our exact 0{as) results with those of a widely used 
effective Lagrangian approach, the A;, approximation, is also presented. 
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In the MSSM, the production mechanisms for the Higgs bosons can be significantly different from that in the 
Standard ModeL For large values of tan/3, the heavier Higgs bosons, A and H, are predominantly produced 
in association with b quarks. Even for tan/3 5, the production rate in association with b quarks is similar to 
that from gluon fusion for A and H production ( Dittmaier et al.l [2011 |). For the lighter Higgs boson, h, for 
tan /3 > 7 the dominant production mechanism at both the Tevatron and the LHC is production with b quarks 

for light Ma (< 200 GeV) , where the bbh coupling is enhanced . Both the Tevatron ( Benjamin et al.l |2010l |) 
and the LHC experiments (jChatrchvan et aP |201l| ) have presented limits Higgs production in association with 
b quarks, searching for the decays h —>■ t'^t~ and bb. These limits are obtained in the context of the MSSM are 
sensitive to the 6-squark and gluino loop corrections which we consider here. 

The rates for hh associated production at the LHC and the Tevatron have b e en extensively studied 

ti 

m the 4- and 5- flavor number parton schemes 
In the 4- flavor number sc heme, the low- 



Dawson et ah" [2006^ , 'Campbell et al.' '2004"| , 'Maltoni et al 
Dicus et al. [1999], Da wson et al. 2004 ], Maltoni et al. 2005, 
and the NLO QCD correction are well un d erstood, both 
(iDawson et al.l |2006| . ICampbell et al.l |2004l| . iMaltoni et al.l [20031 



'20031, 'Dawson et al 
, Campbell et al. 



^ '20051 .'Dittmaier et al.l 12004 
2003], Carcna et al. 20071 fl999 



est omer processes tor produ c mg a tiiggs poson ana a quaru are gg — > bbh and qq — >■ oon i|tJawson et al.l 
[2005j l. iDittmaier et al.l [2004 1. iDicus et al.l [l999| ). In the 5— flavor number scheme, the lowest order process 
is hg — >■ hh (hg — >■ hh). The two schemes represent different orderings of perturbation t heory and calculations 
in the two schemes produce rates which are in qualitative agreement ( Dittmaier et al.l [2011 [. ICampbell et al.l 
200|). In this paper, we use the 5-flavor number sche me for simplicity. The re s ummation of threshol d log- 



hhh (Dawson et al.l 



arithms (|Field et al.l I2007II'). electroweak corrections ((Dawson and Jaiswall [2010j l. iBeccaria et al.l [201(]| |) and 
SUSY QCD corrections ( Dawson and JacksonI [20081 ]]) have also been computed for bh production in the 5- 
flavor number scheme. 

Here, we focus on the role of squark and glui no loops. The properties of the S U SY QCD correction s 
to the bbh v e rtex, both for the decay h — 



Guasch et al.' f2003 l) and t he production, bb 



bh (iDabelsteinI Il995 l 
h (Dittmaier et al 



Hall et al.l 119941. ICarena et al.l l2000l 
2004| . iGuasch et al.l [2003j . iHaber et'all 



^001], Harlander a nd Kilgo rc 2003]), were computed long ag o. The contributions from h squarks and gluin os 
to the lightest MSSM Higgs boson mass are known at 2-loops ( Heinemever et all [2005 ]. iBrignole et al.l [20021 ] ]). 



while the 2-loop SQCD contributions to the bbh vertex is known in the li mit in which the Higgs mass is much 
smaller than the squark and gluino masses ( Noth and Spiral [2010l l2008l ]). The contributions of squarks and 
gluinos to the on-shell hbh vertex are non- decoupling for heavy squark and gluino masses and decoupling is only 

achieved when the pseudoscalar mass, Ma^ also becomes large^ 

An effective Lagrangian approach, the Af, approximation ( Hall et al] 19941 ] . ICarena et al.l [2000j ]). can be 
used to approximate the SQCD contributions to the on-shell bbh vertex and to resum the (a^ tan/^/Msf/gy)" 
enhanced terms. The numerical accuracy of the Ab effective Lagrangian approach has been examined for a 
nu mber of cases. The 2— loo p co ntributions to the ligh test MSSM Higgs boson mass of 0{ai,as) were computed 
bv lHeinemever et al.l [20051 ] and IBrignole et al.l [20021 ] . and it was found that the majority of these corrections 
could be absorbed into a 1— loop contribution by defining an effective b quark mass using the A(, approach. 
The sub- leading contributions to the Higgs boson mass (those not absorbed into A;,) are then of C(l GeV). 
The A j) approach als o yields an excellent approximation to the SQCD corrections for the decay process h ^ bb 



( Guasch et al.l 2003I ]). It is particularly interesting to study the accuracy of the Af, approximation for production 
processes where one of the h quarks is off-shell. The SQCD contributions from squarks and gluinos to the 
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inclusive Higgs produ ction rate in association with 6 quarks has been studied extensively in the 4FNS by 
iDittmaier et alj [2007 1. where the the lowest order contribution is gg ^ bbh. In the 4FNS, the inclusive 
cross section including the exact 1-loop SQCD corrections is reproduced to within a few percent using the Af, 
approximation. However, the accuracy of the A;, approximation for the MSSM ne utral Higgs boson production 
in the 5FNS has been studied for only a small set of MSSM parameters in Ref. iDawson and JacksonI [2008j . 
The major new result of this paper is a detailed study of the accuracy of the Af, approach in the 5FNS for the 
bg — > bh production process. In this case, one of the b quarks is off-shell and there are contributions which are 
not contained in the effective Lagrangian approach. 

In this article, we give a brief review of the effecti ve Lagrangian approx imation in section 1. In section 2, 
we summarize the SQ CD calculation s for b g — >■ bh (jPawson and JacksonI [2008]) including terms which are 
enhanced by nib tan f3 ( Dawson et al.l 201lj ). Analytic results for the SQCD correc t ions to bg — >■ bh in the 
extreme mixing scenarios in the b squark sector have been calculated bv lDawson et all |2011 | and are presented 
in Section 3. Section 4 contains numerical results for the y/s — 7 TeV LHC. Finally, our conclusions are 
summarized in Section 5. 



2. SQCD Contributions to gb bh 

2.1. Ab Approximation: Tlie Effective Lagrangian Approach 

Loop cor r ection s which are enhanced by powers of as tan /3 can be included in an effective Lagrangian approach 
(iHall et al.l |l994j . ICarena et al.l |2000l | . ICuasch et al.l |2003j V Using the effective Lagrangian, which we term 



the Improved Born Approximation (or A^ approximation), the cross section is written in terms of the effective 
coupling, 



= ^^''^'^ ( TTa; ) tanjtana j ' 



where 



, sina\ mbj^iB.) 

gtbh = - ^ [2) 

\C0S(3J VSM 



and the 1-loop contribution to A^ from sbottom/gluino loops is ( Hall et all [1994 1 . ICarena et"aL I |2000lll994j ') 



Afc = -^f^Mg^i tan /3/(Af,- , M,~^ , Mg) , (3) 

where the function /(a, 6, c) is. 

The Improved Born Approximation consists of rescaling the tree level cross section, ctq, by the coupling of Eq. 

m 

™=f^V^o. (5) 
\gbbh/ 

The Improv ed Born Approxim a tion has been shown to accurately reproduce the full SQCD calculation of 
pp — ^ tbH~^ ( Berger et al.l |2005j | . lOittmaier et al.l |2009l |). The one- loop result including the SQCD corrections 



for bg — >■ bh can be written as, 

l^SQCD = CT/BA ( 1 + AsQCD ) , (6) 



where Asqcd is found from the exact SQCD calculation summarized in Appendix B of ( Dawson et al.l [201 1| ). 
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Figure 1; Feynman diagrams for g{qi) + b{q2) — )■ b(pi,) + h{ph)- 



2.2. Full One-loop SQCD Contributions to gb bh 



The SQCD contributions to the gb — b h process have been computed in Ref. ( Dawson and JacksonI |2008j ) 
in the mi, — hmit and further, in Ref. ( Dawson et ahl |2011 |) where terms which are enhanced by mf,tan/3 
have been included. 

The tree level diagrams for g{qi) + 6(92) ^ b{Pb) + ^iph) are shown in Fig. [T] The amplitude can be written 
as a sum of following dimensionless spinor products 



Aff = q!^ 



u{pb)u{q2) 



u{pb) I'^u (92) 
mi, 

f, u{pb)iiu{q2) 

mbt 
iu{p b)iiu{q2) 
mbS 



= 92- 



(7) 



where s = (qi + 92)^, ^ = {Pb ^ li)^ and u = [pb — 92)^- In the mf, — limit, the tree level amplitude depends 
only on and M^^, and is generated at one- loop. When the effects of the b mass are included, M^*, , 



and M4 are also generated. 
The tree level amplitude is 

-4o/3lo = -gs9bbh{T'')^pe^{q,){Mi: + M^'} 
and the one loop contribution can be written as 



47r 



-Qsgbbh 



(8) 



(9) 



For detailed calculation of counter-terms and the coefficients Xj, cf. ( Dawson et al.l [201 



3. Results for Maximal and Minimal Mixing in the &-Squark Sector 
3.1. Maximal Mixing 

The SQCD contributions to bg — > bh can be examined analytically in several scenarios. In the maximal 
mixing scenario, 

\^l-rhl\«j^^\Xi,\ . (10) 
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We expand in powers of 



Mt 



In this case the sbottom masses are nearly degenerate, 
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This scenario is termed maximal mixing since 
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(TO|-TO|^/(1 + Ab)^ 



(11) 



(12) 



We expand the contributions of the exact one-loop SQCD calculation (see Appendix B of iDawson et al.l |2011 |) 

Ml 



in powers of l/Mg, keeping terms to O 



Mi 



and assuming Ms 



Ms 



Ah 



nriR » 



Mwt Mz, Mh Mew- In the expansions, we assume the large tan/? limit and take TObtan/? ~ 0{Mew)- 
This expansion has been studied in detail for t he decay h — > bb, w ith particular emphasis on the decoupling 
properties of the results as Ms and Mg — )• oo ( Haber et al. l200lh. Re cently, we studied this expansion of 
SQCD contributions to the process, bg — ^ bh in Ref ( Dawson et al.l |2011 [). Here, we present the final results 
of the calculation. The amplitude squared, summing over final state spins and colors and averaging over initial 
state spins and colors, including one-loop SQCD corrections is 



1:41^ 



Q \9bbh) 



27r 




(13) 



where, is the rescaled bottom quark Yukawa coupling and the subleading terms 

of 0{Mj^y^ /M'g) are given in ( Dawson et al.l [201l|). These subleading terms are not included in the IB A 
approximation and their conribution is usually small. 



Sgbbh 

dbbh 



(2) 



ana on^nax 



3.2. Minimal Mixing 

The minimal mixing scenario is characterized by a mass splitting between the b squarks which is of order the 
b squark mass, | M? — M? |~ M|. In this case, 

mf — Tor >> — — , 

and the mixing angle in the b squark sector is close to zero. 



cos 26h ^ 1 



{M? -M?)^\l + Ab 



1 



As in the previous section, the spin and color averaged amplitude-squared is. 
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The contributions which are not contained in o'jba, 



Sgb, 



Qbbh 



Ms 

(2) 



(14) 



(15) 



(16) 



and (5KmmT7T are given in ( Dawson et al.l 



[201 1( ) and again found to be suppressed by O 



Mew 
Ms 
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4. Numerical Results 



The numerical results for pp — >■ 6(6) /i at y/s = 7 TeV were presented in ( Dawson et alJ 201 The renormal 



ization an d factorization scales w ere chosen to be ^ji = = Mh/2 and the CTEQ6m NLO parton distribution 
functions ( Nadolskv et alJ 2008j |) were used. Figs. [U [3] and S] show the percentage deviation of the complete 



one-loop SQCD calculation from the Improved Born Approximation of Eq. [5] for tan /3 = 40 and tan /? = 20 
and representative values of the MSSM parameters. In both extremes of 6 squark mixing, the Improved Born 
Approximation approximation is within a few percent of the complete one-loop SQCD calculation and so is a 
reliable prediction for the rate. This is true for both large and small Ma- In addition, the large Ms expansion 
accurately reproduces the full SQCD one-loop result to within a few percent. These results are expected from 
the expansions of Eqs. 1131 and [T51 since the terms which differ between the Improved Born Approximation and 
the one- loop calculation are suppressed in the large Ms limit. 

Fig. [3] compares the total SQCD rate for maximal and minimal mixing, which bracket the allowed mixing 
possibilities. For large Ms, the effect of the mixing is quite small, while for Ms ~ 800 GeV , the mixing effects 
are at most a few fh. The accuracy of the Improved Born Approximation as a function of mn is shown in Fig. 
|5]for fixed Ma,^., and tol. As niR is increased, the effects become very tiny. Even for light gluino masses, the 
Improved Born Approximation reproduces the exact SQCD result to within a few percent. 



Maximum Mixing Scenario, LHC 7 

lan P=40, n=200 GeV, M^=Mj=l TeV 



Exact SQCD 

Large M^, expansion 




'600 800 1000 1200 1400 



Figure 2: Percentage difference between the Improved Born Approximation and the exact one-loop SQCD calculation of 
pp — >■ bh for maximal mixing in the 6-squark sector at ^/s — 7 TeV, tan /? = 40, and Ma = 1 TeV. 



In Fig. [71 we show the scale dependence for the total rate, including NLO QCD and SQCD corrections 
(dotted lines) for a representative set of MSSM parameters at y/s = 7 TeV. The NLO scale dependence is 
quite small when jiji ~ ^ Mf^. However, there is a roug hly ^ 5% difference between the predictions found 
using the CTEQ6m PDFs and the MSTW2008 NLO PDFs'M artin et all j2009l |. In Fig. H we show the scale 
dependence for small ^p (as preferred bv Maltoni et al. |200 5jl), a n d see that it is significantly l arger than in 
Fig. [71 This is consistent with the results of Dittmaier etld] [2011 1 . iHar lander and Kilgord [2003 1. 



5. Conclusion 

The analytical and numerical results presented in the previous sections clearly demonstrate that deviations 
from the Af, approximation are suppressed by powers of {Adpw /Ms) in the large tan/3 region. The A;, approx- 
imation hence yields an accurate prediction in the 5 flavor number scheme for the cross section for squark and 
gluino masses at the TeV scale. 



References 

S. Dittmaier et al. (LHC Higgs Cross Section Working Group) (2011), 1101.0593. 



6 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 




Figure 3: Percentage difference between the Improved Born Approximation and the exact one-loop SQCD calculation of 
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Figure 4: Percentage difference between the Improved Born Approximation and the exact one-loop SQCD calculation 
for pp — >■ bh for minimal mixing in the b squark sector at yfs = 7 TeV . 
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Figure 6: Percentage difference between the Improved Born Approximation and the exact one-loop SQCD calculation 
for pp bh as a function of mn at ^/s = 7 TeV and tan /3 = 40. 
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Figure 7: Total cross section for pp — > b{b)h production including NLO QCD and SQCD corrections (dotted lines) as 
a function of renormalization/factorization scale using CTEQ6m (black) and MSTW2008 NLO (red) PDFs. We take 
Mg = 1 TeV and the remaining MSSM parameters as in Fig. [2l 
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Figure 8: Total cross section for pp -> h(b)h production including NLO QCD and SQCD corrections as a function of the 
factorization scale using MSTW2008 NLO PDFs. We take Mg = 1 TeV and the remaining MSSM parameters as in Fig. 
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